The effect of cusps in time-dependent quantum mechanics 
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Spatial cusps in initial wavefunctions can lead to non-analytic behavior in time. We suggest 
a method for calculating the short-time behavior in such situations. For these cases, the density 
does not match its Taylor-expansion in time, but the Runge-Gross proof of time-dependent density 
functional theory still holds, as it requires only the potential to be time-analytic. 
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Within the Born-Oppenheimer(BO) approximation, a 
vast amount of useful information (such as geometry, 
thermodynamics, and vibrations) of a molecule or solid 
can be extracted from the electronic ground state. But 
electronic excitations are important for many areas, from 
photochemistry to photocmission spectra[l]. Ground- 
state density functional theory(DFT)[2-4] has been very 
successful for the former problem, and its time-dependent 
analog (TDDFT)[5] has become popular for the latter [6]. 
Within linear response, TDDFT yields useful predictions 
for the excited states of many molecules [1], and exten- 
sions to solids are a keen area of research [7, 8]. 

All density functional theories rely on a one-to-one 
correspondence between the one-body potential, such as 
— Z/r for an atom, and the density n(r), under some 
restrictions. For ground-state DFT, the particle statis- 
tics and interaction are fixed. For the time-dependent 
case, one also specifies the initial wavefunction. The 
original proof in the ground-state case of Hohcnbcrg and 
Kohn[2] has been refined over the decades [9-11], but its 
essence remains unchanged. In the time-dependent prob- 
lem, after pioneering works by others[12-18], Runge and 
Gross [5] gave a proof assuming the one-body potential 
is time-analytic, i.e., equals its Taylor expansion in time 
around the initial time, for a finite time-interval. Despite 
recent attempts[19-23], no generally applicable proof has 
been found that avoids this expansion. 

Modern DFT calculations employ the Kohn-Sham 
scheme [3], in which a fictitious set of non- interacting elec- 
trons reproduces the one-electron density of the real sys- 
tem. If we assume that time-analytic potentials yield 
time-analytic densities, van Lccuwcn[24] gave a construc- 
tive procedure for finding the TD KS potential for a 
given system, a problem that has not yet been gener- 
ally solved for the ground-state case. But, in the usual 
treatment, matter has cusps in its ground-state electronic 
wavefunctions at the nuclei, and the resulting spatial non- 
analyticities are coupled to the time-dependence in the 



TD Schrodinger equation[21]. Even in the most mundane 
example, a hydrogen atom in a suddenly switched elec- 
tric field, this coupling leads to non-analytic short time 
dynamics, and its density is not time-analytic. 

We develop a method for extracting the exact non- 
analytic short-time behavior of the Schrodinger equation 
in the presence of cusps. There are distinct spatial re- 
gions with different asymptotic behavior for short times, 
and by 'asymptotic behavior' we mean a series expan- 
sion around the initial time t = 0, whose error vanishes 
as t — > 0+. We calculate the exact short-time behav- 
ior for a hydrogen atom in an electric field, and demon- 
strate agreement with linear response theory in the limit 
of weak fields. The constructive procedure for the TD 
KS potential fails for our examples. (In DFT-spcak, the 
u-rcprcsentability question is not solved.) Nonetheless, 
the original proof of Runge and Gross remains valid. 

Consider a single particle prepared in a field-free 
ground-state, ip , placed in a static field which is turned 
on at t=0, and remains on indefinitely. The exact TD 
wavefunction is 



(1) 



where ej and 4>j(r) are the eigenvalues and functions of 
the Hamitonian operator H, and Cj = (<f>j\tpo). We define 
the time- Taylor series 

^ TE (r,t) =J2c P (r)t p , c p (v) =Y^c n ^^-^ n {T), 

p=0 n 

(2) 

which is the result of the usual practice of Taylor- 
expanding (TE) the time-evolution operator: 



U(t) = l-iHt- H l t z /2 



(3) 



In many cases, the solutions agree, but not necessarily 
when the wavefunction has non-analyticities in space. Al- 
though -0 TE formally solves the TDSE, it may not be a 
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valid solution. Holstein and Swift [25] reported failure of 
the TE solution in a model system, in which tpo has com- 
pact support (all space-derivatives of i/jq vanish at the 
boundary). However, such cases do not occur in routine 
electronic structure calculations, and may not be of con- 
cern in practice. For the rest of this paper, we focus on 
the ubiquitious cusp in ipo{ r ) at a nucleus. 

Begin with a simple example. Start from the ground- 
state wavefunction of the hydrogen atom, 



£3/2 

^o(r) = —=exp(-Zr). 



(4) 



Our potential for t > is simply zero: The nucleus has 
been instantaneously vaporized. The exact TD wave- 
function can be found by applying the free-particle TD 
Green's function to ^(0), yielding 

Z 3/2 iZ 2 t/2 

^(r,t) = —-= [f(r,t) - f(-r,t)}, (5) 

where f(r,t) = (r + iZt) exp(Zr) erfc[(r + iZt)/y/2i£\. 
Fig. 1 shows the exact solution and the power-series 
solution as in Eq. (2), which can be summed to all orders: 

7 3 / 2 , 7t 

V n T 

The TE density diverges as r — > and so is invalid. 
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FIG. 1. TD exact density and TE density of the ground-state 
wavefunction of hydrogen atom under free-propagation. 

Our method begins with introducing reduced variables: 
Ze 2 



m e t 



Aiteoh V h ' 




(7) 



where the powers are motivated by the previous exam- 
ple, and constants chosen for convenience. These reduced 
variables are appropriate if the system contains only one 



nucleus placed at the origin. In reduced variables and 
atomic units (e 2 = h = m = 1), the TDSE becomes 



where 



{£ - isd/ds + 2s 2 v/Z 2 } ip = 0, 



C = -V 2 /2 + zf • V, 



(8) 



(9) 



Z is the nuclear charge, and v is the external potential. 
In the disappearing-nucleus case, ip is equal to its Taylor- 
expansion in powers of s for fixed f , and thus we try: 



p=0 



(10) 



yielding: 



{£ - mi}ip m (r) 



^u p _ 2 (r)V>m- P (r), (11) 



p=0 



where v(r, i) = J2p v p(*) sP ■ Each half-power of t pro- 
duces a second-order differential equation for a function 
of r. The initial wavefunction in reduced variables is 
i/jo(y/2fs/Z). For finite values of r, as s — > 0, f — > oo, so 
its expansion as s — > determines the large f behavior 
of the Vm(f). 

To find the leading non-analytic behavior in t due to a 
cusp, we rearrange the TE solution of Eq. (2) in terms of 
the reduced variables s and f , and test its validity power 
by power in s. Each '0™ must satisfy the differential 
equation, the boundary condition derived from the ini- 
tial wavefunction as f — ¥ oo, and remain finite at the 
nucleus (m being the order of the s-expansion) . We do 
this recursively, until we find k, the lowest value of m for 
which V'fc E fails. Define 



C(f, S ) = ^(r, S )-V it (f ) s) 



(12) 



and solve the differential equation for £fc(r) to find the 
leading non-analytic behavior due to the cusp. 

This can be quite demanding as, in 3D, the equation 
is a partial differential equation. A simple approach is to 
perform a local analysis using the method of dominant 
balance [26] to generate an asymptotic series in inverse 
powers of f . Repeated application yields the series to all 
orders, but this is insufficient to ensure that ip remains fi- 
nite at the r = 0. This difficulty can be overcome by per- 
forming a Borcl resummation of the asymptotic series [26] 
to find the exact solution. The conditions where Borel 
resummation works are discussed in [27]. 

For the hydrogen atom in a suddenly switched electric 
field which stays on indefinitely, the initial wavefunction 
is that of Eq. 4, and the potential after t = is 



v(r) = —Z/r + £ , 



(13) 



where £ is the amplitude of the electric field. This situa- 
tion is covered by the Runge-Gross theorem as a potential 
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that is analytic in time (in fact, constant) with an initial 
wavcfunction that is not its instantaneous ground state. 
The time-dependent dipole moment /i(t) determines the 
photoabsorption of the atom via 



/ dr 
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dt 


t=T- 



sin(o;r), 



(14) 



known analytically[28, 29] for the H atom as £ — >• 0. 

To apply our method, we first calculate the TE wave- 
function order-by-ordcr. We find each term satisfies both 
the differential equation and boundary conditions until 
the 4th order: 



^JE = ^ =( _ 3 + 12z -2 +4f 4 ) 



240F V 



£ 



12^Z 3 / 2 i 



■(l + 6if 2 + 24f 4 ) (15) 



which works except that it diverges at the nucleus. As- 
sume a correction of the form: 



£ 4 = exp(S(r))z, 



(16) 



because the applied field is linear in z. This yields an 
ordinary differential equation for S: 



S" + (S") 2 - 2ifS' + AS' If + 6i = 0, 



(17) 



and applying the method of dominant balance recur- 
sively, one finds the asymptotic expansion: 



ci f" + 



9if 



9 

4f 



3i 
8r* 



C2- 



1 



^ oo 
i2 X] 



(-i) m+1 (m + 4)(2m + 6)! 
(m + l)!2 2m + 5 • 9f 2m 



(18) 



Since 0j E satisfies the boundary condition at large f , £4 
must vanish here, and c\ = 0. But this asymptotic series 
does not yield the small f limit. To find this, we perform 
a Borcl rcsummation, yielding 



C2 



1 + t 

72f 3 



(2 + 2i)e lf f (-3 + l&if 2 + Af 4 ) 



2n(3i - 18f 2 + mr + 8f b ) erfc 

Requiring 04 to remain finite yields 

c 2 = (1 - i)£/(2nZV 2 ). 



V2 



(19) 



(20) 



One can easily check analytically that 04 = 0j E + £4 is 
then the unique solution to Eq. (11) for k = 4. 

With -04, we find the leading half-power of the wave- 
function by changing the variables back to (r, t) and tak- 
ing the asymptotic expansion around the t = 0. The 
leading time-half-power is then 



ii)£Z 



5/2 p ir 2 /2t 



■t 



11/2 



> Vi, (21) 



and one can show 4} m >4 contributes to higher order time- 
half-powers by the method of dominant balance. This 
is also the lowest order in t at which rapid non-time- 
analytic oscillations of the phase appear. 

Our derivation applies for any value of £. The one- 
body potential appears linearly in our equations for 
n (r), and only for n > 3. Thus if 04(f) produces the 
leading non-analytic behavior in t, then this behavior is 
linear in £, and is determined exactly by linear response 
theory. (Higher-order effects in £ produce changes in 
other terms, but not the leading non-analytic short-time 
behavior.) To see this, take the Lehmann representation 
of the Green's function of the hydrogen atom and evalu- 
ate the change in the wavefunction to first order in £ . We 
did not find a closed-form, but the stationary phase ap- 
proximation to the integral over the wave number yields 
the leading half-power term in agreement with our short- 
time result Eq. (21). 

The leading half-power term in the induced dipole mo- 
ment calculated using linear response theory yields: 



fj,(t-> 0+) 256Z 5 /(2835v^) t 9/2 



(22) 



Its Fourier transform yields the known high-frequency 
decay of the photoabsorption cross-section[28] 



a(u — s- 00 ) - 16V2Z 5 tt/(3 cu> 1 . 



(23) 



Thus the cusp at the nucleus leads to fractional powers in 
time-evolution and fractional powers in frequency decay. 

The leading half-power term in a(to — > 00) calculated 
with 04 differs from Eq. (23) by a factor of 2. In calcu- 
lating a(uj) from ip(r,t), one must integrate over r, and 
the coupling between r and t in (as seen in 04 ) allows 
higher order terms in the s-expansion to contribute to 
the leading half-power term in a(uj — > 00) as well. In this 
case, 04 and 05 both contribute to the oj~ 7 / 2 term in a. 

Before discussing the consequences for TDDFT, we 
again consider the simpler, but starker example of the 
disappearing nucleus case. Our method correctly shows 
the error appearing at s 1 order, and produces the correct 
short-time behavior: 

r . Z*l 2 er Zr iZ^ 2 (Zr-2)e- Zr 
0(|r| > Vt., i -> 0+) = — + V , ^ ; 1 



Z 9 / 2 {Zr~4)e- Zr t2 _ (2 + 2i)Z 5 /^ 2 /^ 5/2 
8y/nr irr 4 



(24) 



while its TE contains only integer powers of t. 

Although less relevant for TDDFT, our method applies 
equally well in Id and proves that the short-time behavior 
for the initial state with cusp example in [21] originates 
from s 1 order and has fractional power i 3 / 2 , again with 
a rapid non-time-analytic phase. Similarly, it's s 4 order 
for a particle in a delta-well with an applied electric field, 
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and the leading short-time behavior in ip is 

_ (4 + ^)^/^^/(20]^ ^ 

consistent with the known polarizability and its high- 
frequency limit. [30] 

Our results are proven only for single electrons. How- 
ever, the nuclear potential dominates over the electron- 
electron repulsion in the region of the nucleus, so that 
the cusp condition on the time-dependent density near a 
nucleus remains valid regardless of the interelectron re- 
pulsion. Thus we expect the qualitative features (i.e., 
half-powers of t in the short-time behavior) to remain 
true even in the presence of interaction. 

We conclude by discussing the relevance of our results 
for TDDFT, whose theorems certainly apply (or not) for 
N = 1. The proof of a one-to-one correspondence [5] 
shows that, for two t- analytic potentials whose Taylor 
expansions first differ in the fc-th order, the difference in 
the k + 2th-order t-derivative of the densities is non-zero. 
Thus the two densities differ, the Rungc-Gross proof ap- 
plies, and the potential is a functional of the density un- 
der the conditions stated [5]. Whether or not the TE 
density matches the true density is irrelevant. 

On the other hand, the KS potential of Rcf. [24] pro- 
duces a density whose i-derivatives equal those of the 
interacting density. If we assume that the density is v- 
representable by a i-analytic potential (e.g., if the KS 
and interacting wavefunctions themselves are i-analytic), 
then that potential will yield this density. But the wave- 
functions are not i-analytic in the examples given here, 
and typically are not when an initial wavefunction with 
a cusp undergoes a non-trivial evolution. In general, 
a more sophisticated procedure is needed to ensure the 
constructed potential generates the desired density, and 
some progress has been made [21-23]. This is important 
as the constructive procedure has been invoked or applied 
as is to a variety of situations[31-35]. 

One might argue that real atoms have finite nuclei, 
or that real molecules and solids have nuclear wavefunc- 
tions that smear the cusps due to nuclei. But such ar- 
guments miss the basic point. Just as in ground-state 
DFT, time-dependent DFT is an exact mapping of the 
quantum mechanics of electrons, for which there are no 
difficulties within BO or need for finite nuclei. 
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